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Nonlinear  Singular  Sturm-Liouville  Problems 
and  an  Application  to  Transonic  Flow  Through  a  Nozzle 


SZE-BI  HSU 


Institute  of  Applied  Mathematics  National  Tsing  Hua  University 


We  consider  a  class  of  singular  Sturm-Liouville  problems  with  a  nonlinear  convection  and  a 
strongly  coupling  source.  Our  investigation  is  motivated  by,  and  then  applied  to,  the  study  of 
transonic  gas  flow  through  a  nozzle.  We  are  interested  in  such  solution  properties  as  the  exact  number 
of  solutions,  the  location  and  shape  of  boundary  and  interior  layers,  and  nonlinear  stability  and 
instability  of  solutions  when  regarded  as  stationary  solutions  of  the  corresponding  convective 
reaction-diffusion  equations.  Novel  elements  in  our  theory  include  a  priori  estimate  for  qualitative 
behavior  of  general  solutions,  a  new  class  of  boundary  layers  for  expansion  waves,  and  a  local 
uniqueness  analysis  for  transonic  solutions  with  interior  and  boundary  layers. 


1.  Introduction 

Consider  the  nonlinear  Sturm-Liouville  problem 

eu"  =  f(x,  u)'  -  c{x)h(u ) 

with  u  prescribed  at  x  =  0  and  x  =  1.  We  study  properties  of  the  solutions:  their 
number,  their  asymptotic  shape  for  small  e,  and  their  stability  and  instability, 
when  viewed  as  stationary  solutions  of  the  corresponding  time  evolution  equa¬ 
tion.  Our  study  is  motivated  by,  and  then  applied  to,  the  problem  of  transonic 
gas  flow  through  a  nozzle.  Novel  elements  in  our  theory  include  a  new  class  of 
boundary  layers  corresponding  to  expansion  waves,  a  priori  estimate  for  qualita¬ 
tive  behavior  of  general  solutions,  and  a  local  uniqueness  analysis  for  transonic 
solutions  with  interior  and  boundary  layers. 

For  simplicity  of  presentation,  we  shall  carry  out  our  analysis  for  the  model 
problem 

(1.1)  eu"  =  f(u)'  -  c(x)h(u),  u  =  u(x)  G  R1.  Osisl, 

(1-2)  u(0)  =  u  i ,  w(l)  =  ur. 
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The  flux  function  /(«)  is  assumed  to  be  convex.  This  is  motivated  by  gas 
dynamics,  where  the  sound  speed  depends  monotonically  on  the  density.  By 
composition  with  a  simple  translation,  we  may  assume  that 


(1.3) 


f"(u)  >  0  for  all  u  under  consideration, 
/(?)  =/'(0)  =0. 


The  function  li(u)  represents  the  coupling  of  the  source  due  to  the  geometry  and 
the  gas  flow.  The.  following  strong  nonlinear  coupling  assumption  is  dictated  by 

physics:- 


(1.4)  h(u)  *  0.  h'(u)  0  for  all  u  under  consideration. 


The  function  <  (  )  represents  the  strength  of  the  source  and  may  change  sign.  For 
simplicity  and  definiteness,  we  make  the  following  assumption: 

(1.5)  i(v)/t(u)<0  for  0  g  v  <  1  and  all  u  under  consideration. 

More  general  situations  can  also  be  treated  with  our  analysis.  There  are  two 
distinct  cases: 

( 1 .6)  i  (  v  )li'(  u  )  <  0  (stable  case;  diverging  duct). 

(1.7)  <  (  v  )li'(  u )  >  0  (unstable  case;  converging  duct). 

Our  equation  (1.1)  may  he  regarded  as  the  stationary  equation  for  the 
convective  reaction-diffusion  equation 

(1.8)  u,  +  f(u)  x  =  euxx  +  c(x)h(u). 

In  the  stable  case.  (1.1)  (1.2)  has  a  unique  solution  which  is  time  asymptotically 
stable  with  respect  to  (1.8).  In  the  unstable  case,  there  may  exist  three  solutions, 
one  of  which  has  an  interior  layer  and  is  nonlinearly  unstable  with  respect  to 

(1.8) ;  the  other  two  are  stable. 

We  now  relate  our  problem  to  nozzle  flow.  The  quasi-one-dimensional  model 
for  isentropic  viscous  compressible  flow  through  a  nozzle  is 

<  1  -*,)i  ipA),  +  (prA)x  =--  0. 


( 1  .^):  (  prA),  4  (pr:A  )  x  +  Apx  =  p(  Ai\  )  v. 

where  p.r.p.p  and  A  =/f(.v)  are  the  density,  velocity,  pressure,  viscosity 
coefficient  of  the  gas  and  the  cross  section  of  the  nozzle.  The  pressure  is  a  given 
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function  of  the  density  p.  We  assume  that 

(1.10)  p'(p)  >  0,  p"(p)>  0. 

This  is  satisfied  by  the  polytropic  gases  where  p(p)  =  apr.  f  §  r  >  1.  The  station¬ 
ary  equations  are 


(pvA)x  =  0, 

(pv2A)x  +  Apx  =  p(Avx)x 


The  first  equation  can  be  integrated  to  yield  a  relation  between  v  and  p.  This  can 
be  used  to  eliminate  one  of  the  dependent  variables  and  so  yield  a  scalar  equation 
for  one  dependent  variable.  If  p  is  eliminated,  we  get 


(1.11) 

(1.12), 

(1.12); 

(1.12), 


(g(c,  x))x  =  fU'„  +  c(x)k(v ), 


.?(«••  x) 


-  M 


/!'(■*• ) 

A( x) ' ' 


'■'(  v ) 


A'(x) 

(  -4(  a  )) 


;  PoPpA  () 


A'(x)  | 
A(x)  I  ■ 


k(r)  =  v. 


where  p„  =  p(0).  r„  =  r(0)  and  A„  =  .4(0).  From  (1.10)  we  have 


X,  (  r.  x  I 

v) 


Pii|’ii.'1»  _  PqEii^q  ,1  PqIqA,,  \  ^  4  ( .v  ) 

A  (  v  )  r2A(x)P\rA(x)  I  ~  /*,4(.v)  ’ 

2Pd,'<i'4(i  ,(  Pni'i)A()  )  |  p„(’()/l(l  |  /  plli\)A0  ^ 

eJo!'’  1 72TTT I  +  ( P3I7T I  M  ~  I  ’  0 


The  first  equation  shows  that,  for  small  viscosity  p.  g,  >  0  if  the  flow  is 
supersonic  (i.e.,  when  |r|  >  ( p')1  :).  while  g,.  <  0  if  the  flow  is  subsonic  (Le¬ 
w-hen  |r|  <  (  p')1  :).  We  are  interested  in  transonic  flow.  i.e..  when  |r|  is  close  to  ,‘ 
the  sonic  speed  (  p')‘- :.  Thus  we  assume  that  (r|  *  0.  and  so  k(r)  *  0.  We  see 
that  the  nozzle  equation  (1.11)  shares  properties  analogous  to  (1.3).  (1.4)  for  the 
model  (1.1).  For  this  reason,  we  will  call  the  positive  (negative)  states  u  >  0 
(w  <  0)  for  (1.1)  supersonic  (subsonic)  and  call  the  zero  state  sonic.  It  is  clear 
from  the  above  equations  that  for  a  converging  duct.  i.e..  rA'(  v )  <  0,  we  have  n~ 
<•(  v ) k '( r)  >  0.  This  is  the  reason  we  called  the  unstable  case  (1.7)  the  converging  ~ 

duct  and  the  stable  case  (1.6)  the  diverging  duct.  - 

_Distributlon/ 
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Our  main  interest  is  in  the  properties  of  solutions  for  small  e.  The  inviscid 
theory,  s  =  0+,  has  been  worked  out  by  Liu  in  [4],  This  is  reviewed  in  Section  2. 
We  also  remark  that  a  new  type  of  boundary  layer  with  algebraic  decay  exists. 
They  correspond  to  rarefaction  waves  connecting  to  the  sonic  state.  Remark  2.1. 
In  Section  3  we  present  an  analysis  which  shows  that  any  general  solution  of 

(1.1),  (1.2)  is  close  to  one  of  the  corresponding  inviscid  solutions.  Our  analysis 
differs  markedly  from  the  usual  asymptotic  analysis  where  solutions  of  viscous 
equations  are  constructed  based  on  the  inviscid  solutions.  In  the  stable  case  (1,6), 
it  follows  easily  from  the  maximum  principle  that  there  exists  at  most  one 
solution.  The  inviscid  theory  predicts  that  in  the  unstable  case  there  are  one  or 
three  solutions,  depending  on  the  data  (1.2).  With  the  a  priori  qualitative 
understanding  obtained  in  Section  3,  we  show  in  Section  4  that  (1.1),  (1.2)  with 
(1.7)  has  exactly  one  or  three  solutions  by  proving  a  local  uniqueness  theorem. 
For  solutions  with  no  interior  layer,  this  is  done  by  generalizing  and  lefining  the 
classical  argument  of  Coddington  and  Levinson  [2].  For  solutions  with  interior 
layer,  a  new  argument  is  introduced  for  the  local  uniqueness  theorem. 

In  Section  5  we  present  a  stability  analysis  for  solutions  of  (1.1)  viewed  as 
stationary  solutions  of  (1.8). 

There  have  been  studies  on  singular  nonlinear  Sturm-Liouville  problems:  see 
[3]  and  references  therein.  However,  these  studies  do  not  consider  models  with 
the  strong  coupling  property  (1.4).  As  we  show  here,  this  coupling  property  has 
the  regularizing  effect  that  viscous  solutions  are  close  to  the  inviscid  solution. 
In  the  absence  of  (1.4),  the  inviscid  theory  often  offers  too  many  solutions,  most 
of  which  do  not  correspond  to  viscous  solutions.  With  (1.4)  there  arise  new 
analytical  difficulties,  some  of  which  are  resolved  here.  It  would  be  interesting  to 
study  more  general  convective  reaction-diffusion  equations  such  as  the  full 
quasi-one-dimensional  nozzle  flow  equations  (cf.  references  in  [4]). 


2.  Inviscid  Theory 

In  this  section  we  review  the  time-asymptotic  states  for 

(2.1)  u,  +  f(u)x  =  c(x)h{u). 

We  then  present  a  brief  account  of  that  theory.  Except  for  boundary  or  interior 
layers,  the  solutions  of  (1.1)  should  satisfy 

(2.2)  f(u)x  =  c(x)h(u) 

as  e  -*  0 , .  Interior  layers  would  tend  to  stationary  shock  waves  ( u  .  u  + )  of  (2.1) 
satisfying  the  jump  condition  and  entropy  condition 

(2.3)  /(u  +  )  =/(«  ).  u+<  u  . 

An  inviscid  boundary  layer  (uhu0)  at  x  =  0  would  correspond  to  stationary 
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waves  of(l.l)atx  =  0  which  in  turn  satisfy  approximately 

f(u)x  =  euxx ,  u(0)  =  tf(oo)  -  u0,  ux( oo)  =  0, 

since  in  the  layer  the  source  c(x)h(u)  has  little  effect.  Integrate  the  above 
equation  from  x  =  x  to  x  =  oo  to  get 

(2.4) 0  eux(x)  =  f(u(x))  -  f(u0),  u(0)  =  uh  «(ao)  =  u0. 

In  order  for  the  above  equation  to  have  a  solution,  we  need 

(2.4)  ( u  —  uQ )( /( u)  -  f(u0 ))  <  0  for  all  u  between  u,  and  u0. 

Thus  an  inviscid  boundary  layer  (uh  u0)  of  (2.1)  at  x  =  0  would  correspond  to  a 
limit  of  viscous  boundary  layers  of  (1.1)  if  and  only  if  (2.4)  holds.  Similarly,  a 
boundary  layer  («0,  ur)  at  x  =  1  satisfies 

(2.5) 0  eux(x)  =  f(u(x)) -/(u0),  u(\)  =  ur.  m(-oo)  =  u„. 

(2.5)  (m -m0)(/(u)- /(u0))  >  0  for  all  u  between  ua  and  ur. 

Conditions  (2.4)  and  (2.5)  can  be  related  directly  to  elementary  waves  for 

(2.6)  «,+/(«)„  =  0. 

There  are  two  types  of  waves  for  (2.6).  Two  states  u  and  u  .  can  be  connected 
by  a  shock  wave  ( u  .  u . )  with  speed 

/(«,)-/(«*  ) 

o  =  — - - - — 

u ,  —  u 


if  u,<  u  .  When  u.>  u  (k_.h  +  )  is  a  rarefaction  wave  propagating  with 
characteristic  speed  f'(u).  This  is  so  under  the  convexity  condition  /"(«)  >  0. 
(1.3).  ft  can  be  seen  easily  that  (2.4)  holds  if  and  only  if  (u,.  w„)  constitutes  an 
elementary  wave  for  (2.6)  with  negative  speed.  Similarly.  (2.5)  holds  if  and  only  if 
(u0.  ur)  has  positive  speed. 

Remark  2.1.  Formal  asymptotic  expansion  for  boundary  layers  at  x  =  1  can 
be  performed  to  yield  (2.5 )„  as  follows:  W  rite 

u(x)  =  U{x)  +  r(  tj). 

1  -  x 


36 


S.-B.  HSU  AND  T.-P.  LIU 


where  L'{  \)  is  valid  outside  the  layer  and  therefore  is  close  to  being  an  inviscid 
solution  satisfying  (2.2),  and  c(t))  is  valid  in  the  layer.  Plug  the  expression  into 
(1.1 )  and  compare  the  coefficients  of  t  to  obtain 

r"  +  f'(c)c'  =  0, 


r(0)  =  ur.  c(x)  =  r'(oc)  =  0. 


Integrate  above  to  yield  (2.5)„. 


r'  =  -fir).  r< 0)  -  ur.  r(oc)  =  0. 


From  this  we  notice  that  in  the  case  of  a  rarefaction  boundary  layer,  ur  >  0.  r ( 77 ) 
decays  algebraically  as  a  consequence  of  (1.3).  When  f{u)  has  higher-order  zeros 
at  u  =  0.  the  decay  rate  is  lower.  This  is  in  contrast  to  the  usual  types  of 
boundary  layers,  where  the  decay  is  exponential. 

The  above  can  be  summarized  as  follows:  Solutions  of  (1.1).  (1.2)  would  tend 
to  inviscid  solutions  which  satisfy  (2.2)  except  for  possible  discontinuities.  The 
discontinuities  at  0  <  x  <  l..x  =  0  and  .x  =  1.  would  satisfy  (2.3).  (2.4)  and 
(2.3).  respectively.  Such  an  inviscid  wave  pattern  is  called  an  asymptotic  state 
because  it  represents  the  large-time  state  of  solutions  with  given  end  states  at 
x  =  ±  x  (see  l.iu  [4]).  We  now  describe  all  the  possible  types  of  asymptotic 
states  with  given  end  states  u,  and  ur. 

Besides  hypothesese  (1.3)  and  (1.5)  we  lurther  assume  that 


(2.7) 


lint 


h(u) 


.  ,  f(u) 


=  0. 


Since  c(  .x  )h{u)  <  0.  (1.5).  a  solution  u(x)  of  (2.2)  moves  toward  the  sonic  slate 
zero  as  x  increases.  Condition  (2.7)  ensures  that  given  a  slate  u  there  always 
exists  a  slate  it  such  that  u  and  u  are  connected  by  a  solution  of  (2.2)  with  values 
5  at  v  =  0  and  u  at  x  =  1 .  Since  (2.2)  is  singular  at  u  —  0.  and  /'(0>  =  0  by 
(1.3),  there  are  two  states  u*  and  u„.  with  «*  >  0  >  each  of  which  is 
connected  to  u  =  0  bv  solutions  of  (2.2).  Given  a  state  u  *  0.  define  u.  uu  <  0. 
satisfying 


/(«)=/(«) 

so  that  u  and  u  form  a  standing  shock  wave,  (2.3). 

Consider  first  the  stable  case  (1.6),  c(x)h'(u)  <  0.  It  is  easily  shown  that  in 
this  case 


(2.8) 


u *  >  u,. 
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C use  A.  u i  g  ut. 

Case  Al.  When  ur  a  0.  the  asymptotic  state  consists  of  a  backward  wave 
(u,.  ut)  at  .v  =  0,  a  subsonic  stationary  wave  connecting  ut  at  x  =  0  and  u  —  0 
at  x  =  1.  and  a  forward  rarefaction  wave  (0,  ur)  at  x  =  1.  The  backward  wave 
(U/.  m»)  is  a  shock  (rarefaction)  wave  when  u,  >  u„(ut<  u,). 

Case  A2.  When  u,  <  0.  the  asymptotic  state  consists  of  a  backward  wave 
(it,.  iir)  and  a  stationary  wave  (ur.  ur). 

Case  B.  Tit  <  ut  <  u*. 

Case  B1 .  When  u,  J>  0.  the  asymptotic  state  consists  of  a  supersonic  station¬ 
ary  wave  (  u i-  ii  )  for  0  g  x  <  ,v0.  a  stationary  shock  wave  (u  .  u  .  )  at  .v  =  x„.  a 
subsonic  stationary  wave  (t/,.0)  for  x0  <  x  <  1  and  a  supersonic  rarefaction 
wave  (0.  t/,)  at  v  =  1.  The  location  x  =  x„  and  the  states  u  ,  u  .  of  the  standing 
shock  wave  are  determined  uniquely  by  the  left  slate  ut. 

Case  B2.  W'hen  ur  <  0  and  ur  <  ti,.  the  asymptotic  state  consists  of  a 
backward  shock  wave  («,.  ur)  and  a  subsonic  stationary  wave  (ur.  ur). 

Case  B3.  When  u,  <■  0  and  ur  a  u,.  the  asymptotic  state  consists  of  a 
supersonic  stationary  wave  (  m;.  u  )  for  0  §  x  <  v„.  a  standing  shock  wave 
(  ii  .  u  .  )  at  v  =  a,,  and  a  subsonic  stationary  wave  (u  ,.ur)  for  x(l  <  v  g  ]. 

Case  C.  u i  §  t/*.  Define  u,  with  u,  =  u,. 

Case  Cl.  When  u,  g  u,.  the  asymptotic  state  consists  of  a  stationary  wave 
(u,.  m,)  and  a  forward  wave  (ux.  ur)  at  x  =  1. 

Case  C 2.  When  ar  <  ux  <  0  and  ur  a  the  asymptotic  state  consists  of  a 
stationary  wave  for  v  *  ,v(l  and  a  stationary  shock  wave  at  ,x  =  xu. 

Case  C3.  When  ur  <  ul  <  0  and  ur  <  uh  the  asymptotic  state  consists  of  a 
backward  shock  wave  (u,,  ur)  at  x  =  0  and  a  stationary  wave  (ur.  ur). 

Next  we  consider  the  unstable  case  (1.7),  c(x)h'(u )  >  0.  In  this  case  we  have 

(2.9)  u,  >  u*  >  0. 


Case  D.  u,  £  u*. 

Case  Dl.  When  u,  }>  0,  the  asymptotic  state  consists  of  a  backward  wave 
( u ut)  at  x  =  0,  a  subsonic  stationary  wave  ( u 0)  for  0  <  x  <  1 .  and  a  forward 
rarefaction  wave  (0,  ur)  at  x  =  1. 
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Case  D2.  When  ur  <  0,  the  asymptotic  state  consists  of  a  backward  wave 
(u,.  ur )  at  x  =  0  and  a  subsonic  stationary  wave  (ur,  ur)  for  0  <  x  g  1. 

Case  E.  u*  <  u;.  Define  ux  by  ul  =  u 

Case  El.  When  ur^0  and  u*  <  u,  <  ut,  there  are  three  asymptotic  states: 

(i)  a  supersonic  stationary  wave  (uh  «,)  for  0  g  x  <  1  and  a  forward  wave 
(u,.  ur)  at  x  =  1; 

(ii)  a  supersonic  stationary  wave  («,,!*_)  for  0  g  x  <  x0,  a  stationary 

shock  wave  (u_,  u+)  at  x  =  x0,  a  subsonic  stationary  wave  (u  +  ,  0)  for 
x0  <  x  <  1  and  a  forward  rarefaction  wave  (0.  ur)  at  x  =  1;  x0,  «  + 

are  uniquely  determined  by  u,; 

(iii)  a  backward  shock  wave  ( uh  u*)  at  x  =  0,  a  subsonic  stationary  wave 
(z/„.0)  for  0  <  x  <  1.  and  a  forward  rarefaction  wave  (0,  ur)  at  x  =  1. 

Case  E2.  When  ur  ;>  0  and  u,  ^  «„  the  asymptotic  state  consists  of  a 
stationary  wave  ( uh  ux )  for  0  g  x  <  1  and  a  forward  wave  (uv  ur)  at  x  =  1. 

Case  E3.  When  ur  <  0,  there  are  three  subcases:  (i)  if  ul  <  ur  <  0,  then 
an  asymptotic  state  consists  of  a  stationary  wave  ( uh  u, )  and  a  forward  shock 
wave;  (ii)  if  u,  >  ur,  then  it  consists  of  a  backward  shock  wave  (u,,  u,)  and  a 
stationary  wave  (ur,  u,);  (iii)  if  ul  <  ur  <  0  and  U ,  >  u„  then  there  also 
exists  an  asymptotic  state  which  consists  of  transonic  stationary  waves  with  a 
stationary  shock  wave  at  x  =  x0,  determined  uniquely  by  u,  and  ur. 

Thus,  given  end  states  u,  and  the  inviscid  theory  yields  three  solutions 
for  Case  El  and  also  (iii)  of  Case  E3.  It  can  be  shown  easily  that  the  above  is 
a  complete  description  of  asymptotic  states  and  that  an  asymptotic  state 
depends  smoothly  on  its  end  states. 


3.  A  Priori  Properties 

As  in  Section  2.  hereafter  besides  (1 .3)— (1 .5)  we  also  assume  that 


(3-1) 


lim  ^ 

l“l— OC  /(“) 


=  0. 


Lemma  3.1.  Any  solution  u(x)  of  (1.1)  belongs  to  one  of  the  following  three 
types : 

Type  I:  u(x)  is  strictly  increasing; 

Type  11:  u(  x )  is  strictly  decreasing; 

Type  III:  u(x)  has  a  unique  critical  point  which  is  an  absolute  minimum. 

Proof:  This  is  an  immediate  consequence  of  the  hypothesis  (1.5), 
c(x)h(u)  <  0. 
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Lemma  3.2.  There  exists  a  positive  constant  M  depending  on  w,  and  ur  and  not 
on  e  such  that  any  solution  of  u(.x)  =  u(x;  e)  of  (1.1),  (1.2)  satisfies 

|m(x)|  <  M,  OsitSl. 

Proof:  The  lemma  holds  trivially  for  monotone  solutions.  Let  u(x)  be  a 
Type  III  solution  of  Lemma  3.1  with  minimum  u(x0)  at  x  =  xQ.  Integrate  (1.1) 
from  x  =  .x0  to  x  =  1.  and  use  (1.5),  ch  <  0.  to  obtain 


(3.2) 


fi/(l)  +f(u(xn))  =/(«r)  -  J\(x)h(u(x))  dx 


=  /(«,)  +  Jl\c(x)h(u(x))\dx. 

xi\ 


Since  h'(u)  *  0  and  u(x)  is  strictly  increasing  for  x0  <  x  <  1.  h(u(x ))  is 
monotone  for  ,x0  <  .x  <  1.  When 

max  h(u(x))  =  h(u(\))  =  h(ur), 

a  *  s  i 


we  have  from  (3.2)  and  u'(l)  >  0  that 


/(“(*«))  <  /(«,)  +  h(ur)  J'c(x)dx 

x 0 


which  is  bounded  independent  of  e  and  so  the  lemma  is  proved.  When 

max  h(u{x))  =  h(  u(x0)). 


we  have  again  from  (3.2)  that 

/(u(x0))  <  f(ur)  +  h(u(x0))  flc(x)dx. 

x0 

This  estimate  and  hypothesis  (3.1)  yield  an  upper  bound  independent  of  e  for 
|u(x0)  |  and  therefore  for  |u(x)  j,  0  g  x  g  1. 

Lemma  3.3.  Suppose  that  u(x),  a  <  x  <  b.  is  a  strictly  increasing  (decreas¬ 
ing)  solution  of  (1.1),  (1.2)  and  that  u(x)  >  -Cel/i  ( u(x )  <  Ce1/4),  a  <  x  <  b. 
for  some  C  >  0.  Then  b  -  a  =  0(l)eI/2  as  e  -»  0+. 
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Proof:  We  consider  only  the  case  when  u(x)  is  increasing;  the  other  case  is 
similar.  Integrate  (1.1)  from  a  to  x  to  obtain 

eu'(x)  -  tu'(u)  =  f(u(x))  -  f(u(a))  -  f  c(  y)h(u(y))  dy. 

J  a 

Since  u  is  increasing,  u{  \)  >  u(a)  for  x  >  a.  and  u(a)  >  —  CV/4.  we  have  from 
(1.3)  that  /( ti(  x ))  -  f(  u(a))  >  -  iCV  ".  Thus  the  above  yields 

tu'(.x)  ;>  -  5fV';  -  f  c(y)h(u(y))  dy  for  a  <  x  <  b. 

Ju 

Integrate  this  from  a  to  b  to  yield 

£u(b)-eu(a)>  -  —  a  )el/1  -  (  f  e(  y)h(u(  y))  dy  d.x . 

Ja 

By  assumption.  c(.x)h(u)  <  0,  and  from  Lemma  3.2.  |u|  <  M.  Thus,  the  above 
estimate  yields 

2  t'M  +  y  C2  f1  2(b  -  a)  s  D(b  -  a )“. 

for  some  positive  constant  D.  The  lemma  follows  immediately  from  this 
inequality. 

In  the  following  two  lemmas  we  study  the  solution  of  (1.1).  (1.2)  outside  the 
layers.  In  all  cases  we  assume  that  x,  -  ,v0  is  of  order  one.  i.e..  x,  -  x„  >  (V  s 
for  sufficiently  large  C.  and  0  g  x„  <  x,  g  1. 

Li  mma  3.4.  Let  u(.x)  be  a  solution  of  (1.1).  (1.2).  Set 

C,  =  3(  max /"(«))[  max  |r(x)A(w)|  . 

'  |M|  J,  A/  M  <>s  <  s  I  ! 

Ms  V 

(f  (  \  )  =  /’(  u(  x  ))  u'(  x  )  -  l  (,v)/i(  m(  a)). 

(i)  Suppose  that  u'(.x)  >  0,  /'(«( a))  <  -r1  4  for  a,,  g  a  g  a,  and 
|<j  (  a, | )  |  <  C‘|f'  4.  Then  |<p(  x )  |  g  C,f'  4.  .v0  g  a  S  a,,  for  e  sufficiently  small. 

(it)  Suppose  that  u'(x)  <  0.  ffu(.x))  >  r1  4  for  .v0  g  A  g  a,  and  |<p(.v,)|< 
C  ,f  '  4.  Then  |<p(  a  )  |  g  Cy1  4.  x„  g  ,x  g  a,,  for  e  sufficiently  small. 

Proof:  We  shall  prove  (i);  (ii)  follows  by  similar  arguments.  Suppose  that 
the  conclusion  fails.  Then  there  exists  x,  x„  g  x  <  a,,  such  that 

|<p(x)|  g  C,el  4  for  a„  <  x  <  a. 
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and  either 

(I) 

<f(x)  =  <r,F1/4. 

<p'(jc)  ?  0, 

or 

(II) 

<p(.v)  =  -CV,/4. 

cp'(  x )  ^  0 

We  treat  (1):  (II)  is  similar.  Relation  (l.l)  is  the  same  as  £«"(*)  =  tp(.x).  With  (I) 
we  have 

M"(.?)  =  C\e  V4. 

/'(m(.v))m”(.v)  +  f"(u(x))(u'(x)) 


-  c'(x)h(u(x))  -  c(x)h'(u(x))u'(x)  ;>  0. 
From  the  hvpothesis  /'( u)  -  r1  4  the  above  yields 
(>  1  :  ^  f"(u(  v  ))(«'( .v))‘ 

(3.3) 


-  < ■'(  V  )/?(  u(  X  ))  -  £•(  .v  )/;'(  u(  x  ))»’(  v  ). 
On  the  other  hand,  from  (I)  we  have 


/'(//(  v  ))m'(.y )  -  c(x)h(u(\ ))  =  (>'  J. 
and  by  Lemma  3.2.  |t/|  <  A/.  Since,  by  hypothesis.  /'(«)  <  -f1  4.  it  follows  that 


M’(  a  ) ;  s  |  niax  i<  (  v  )/i(  n)  |  +  (>'  4  j  r  1  4. 

■  u  ,-_  U 


This  and  (3.3)  yield 


(■,«•  1 

:  <  2  max  f"(  u ) 

max  1  <•(  .v  )/t(  n )  | 

t  ' 

‘  +  (  i 

id  w 

II 1 

1 

l  i«ii  V 

' 

•+  max  |c'(  v )/;(  m)  | 

ii  ■.  «  •_  I 
III'.  V 


+  max  |<  (  \  )h'{  u) )  max  |i  (  \  )/i(  u)  |f  1  4  +  C', 

n  ■.  . ...  i  \  "c  >  i  i 

|nic  W  |u|<  '< 

This  contradicts  the  definition  of  C,  when  e  is  sufficiently  small. 
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Proposition  3.5.  Suppose  that  u(x)  is  a  solution  of  (1.1),  (1-2)  and  is  strictly 
monotone  over  a  subinterval  (x0.  x,)  of  (0,1).  Then  there  exist  C>  0  and  a 
subinterval  (x2,  x3)  of  (x0,  x,)  such  that,  for  sufficiently  small  e. 
M*)|  =  |/'(«(.v))u'(.v)  ~  c(x)/i(«(x))|s  Ct1/4  over  (x2.x3)  and 

|x2  -  x0|  +  |x3  -  xtl  ^  De3/* 

for  some  positive  constant  D  independent  of  e. 

Proof:  In  view  of  Lemmas  3.3  and  3.4,  we  only  need  to  show  that,  for 
strictly  increasing  u(x).  there  exists  x2  with  |<p(x,))g  Cev 4  and  |x3  -  x0|^= 
0(l)e' and,  for  strictly  decreasing  u(x).  there  exists  x3  with  |<p(x'j)]g  Ce 
and  |.\,  -  x,|  =  0(l)ev *  for  any  given  C  >  0.  Consider  the  case  u(x)  is  increas¬ 
ing.  Suppose  that  |  <p ( x )  |  2;  Ce1/4.  x0  <  x  <  x.  for  some  x  >  x„  and  C  >  0. 
There  are  two  cases: 


u) 

T(x)  >  Ce1  4 

for 

X(1  <  X  <  X. 

or 

(in 

qp(x)  <  -  Ce1/4 

for 

1  > 

V 

V 

In 

Case  ( I )  we 

have  from  (1.1),  f«”(x ) 

=  <P(-' 

c).  that  u"(  x 

< 

v  and  so  by 

integration 

u’(  A 

•f 

1 

All 

u'(  x„ 

)  5  Ce  '  4(- 

In 

tegrate  again 

to  obtain 

C  , 

«(x)  -  u ( .v() )  2  yf  ’  4(x- 

Since  u  is  increasing.  2W  2  u(x)  —  u( x0)  and  also  the  above  yields 


x  -  x„ 


=  E. 


Thus  we  have  shown  that  |qr(x2)  |g  Ct1/4  for  some  x2.  xn  <  x:  <  x„  +  E.  This 
proves  the  lemma  for  Case  (I).  Other  cases  are  treated  similarly. 


We  next  studv  boundary  and  interior  layers. 


LKMMA  3.6.  Use  the  same  notations  as  in  Proposition  3.5. 

(i)  When  u\x)  >  0,  xn  <  x  <  x,.  then  u(x)  <  0.  v„  <  x  <  x„  and  x.  can 
be  chosen  so  that  either  tu"{x)  2  C,rl/4  and  jw(x)  -  m(x0)|=  0(  1  )e  for 
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x2  or  eu"(  x )  g  -Qe174  for  x2:  a,  can  be  chosen  so  that 

either  x3  =  Xj  and f'(u(x}))  g  -  C,£1/4  or  x3  <  xL  and  f'(u(x}))  =  -C,e1/4. 

(li)  WTien  w'(x)  <  0,  x0  <  x  <  x,,  then  u(x)  >  0.  x0  <  x  <  x„  and  x2  and 
x3  can  be  chosen  so  that  x2  =  x0  and  either  x3  =  x,  or  x,  <  x;.  In  the  case 
x,  <  Xj  we  have  either  |m(x)  —  «(x,)|^  L>e3/8  for  XjSxS  x,  or  <p(x3)  = 
-  C,i1/4  and /(u(x3))  -  /(w(xj))  S  -  De3/8  /or  some  D  >  0  independent  of  e. 

Proof:  (i)  The  second  part  is  a  direct  consequence  of  (i)  of  Lemma  3.4.  It 
follows  from  Lemma  3.3  that  u(x})  <  0.  Since  u'(x)  >  0  and  u(x)  is  close  to  the 
inviscid  wave,  by  Proposition  3.5  we  have  u(x)  <  -C,  xn  <  x  <  x2.  for  some 
positive  C  independent  of  e.  When  x2  =  x0,  (i)  holds  trivially.  If  x,  >  x0.  then 
by  (i)  of  Lemma  3.4  we  may  assume  that  | <p ( x2 )  |  =  C,el/4  and  |<p(x)  |  s  C,f'  4 
for  xtl  <  x  <  x,.  There  are  two  possibilities: 


(1) 

<f(> 

•t 

O 

All 

for 

x0  <  x  <  x2  and 

<p(x,)  =  C|f!  4 

(II) 

<r(-' 

O 

1 

VII 

for 

v0  <  x  <  x,  and 

T  ( x : )  =  -C,r! 

In  (I).  i/”( x )  >  0.  Since  qr(.v:)  =  C,f1/4  and  u(x2)  <  -  C. 


(x:)h(u(xz))  C>' 

f'(u(x2))  f'(u(x1)) 


for  some  C  independent  of  f.  Thus  u '( x )  <  C'  for  x„  <  v  <  x2.  From  Proposi¬ 
tion  3.5.  |.v;  -  x„|  =  0(l)f3  x  and  so  we  have 

t/(.v:)  -  m(x„)  =  J  'u'(x)  dx  =  0(1  If3  \ 

Case  (II)  corresponds  to  the  case  eu"(x )  g  -C,*7  4  for  x„  <  v  <  x2  in  (i).  This 
proves  (i). 

(ii)  That  v,  -  x0  follows  from  (ii)  of  Lemma  3.4.  If  v,  <  x,.  then  by  (it)  of 
Lemma  3.4  we  may  assume  that  |qp(x3)  |  =  C,f1,4  and  |<p(x)|s  C,f‘  4  for 
\ ,  <  x  <  x,  so  long  as  f’(u{x ))  >  r1  4.  When  q/.v,)  =  C,f'  4  we  can  show,  just 
as  in  Case  I  above,  that  |u(  x )  -  u(x2 )  |  g  Dt'  11  for  x,  g  x  g  v,.  Consider  next 
the  case  qr( x,l  =  -('if1  4  and  <p(x)  g  -C,?1  4  for  v,  g  x  g  x,  and  /'(id  v)) 
•>  r 1  4.  In  the  case  f'(u(x))  g  f1  4  for  x,  irS  x,.  we  see  by  integrating  (1.1) 
that 

ru'(  v,  )  -  m'(  x, )  =  /(  u(.v, ))  -  /(  u(.tj))  +  j  c(x)h(u(x))  dx. 
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From  <f  ( \  j )  =  -C',f'  4  and  /’(u(a,))  i  e1/4  we  have 


4  +  c(x<)h[u(x})) 
/'(«(v,)) 


=  0(  1  )e 


Thus  in  view  of  the  above  estimates  and  Proposition  3.5  we  conclude  that 
/(  m(.V,))  -/(m(.v,  ))  =  c(.v)/t(  u(  a))  dx  -  fw'(-Vi)  +  fu'(.v,) 

>  J  c(x)h(u(x))  dx  +  fu'(  V,) 

=  <9(1  )|  V,  -  .Vjl  +  0(1)f'  4 

=  0(1)  F1  \ 

This  completes  the  proof  of  Lemma  3.6. 

l.i-.MMA  3.7.  Suppose  that  u(x)  is  a  solution  of  (1.1).  (1.2)  and  Inis  a  minimum 
at  v  =  A,,  not  near  the  boundary  x  =  0  or  x  —  1.  i.e..  (  r  *  <  ,v0  <  1  -  Of'  s  for 
sufficiently  large  positive  constant  C.  Then  there  exists  a,.  ,v2.  a,  <  ,v„  <  .v:.  such 
that  <p  t  .v , )  =  -  0|f'  4.  qr(.v2)  =  0,f‘  4.  |.v;  -  a,|  g  Dr'  \  (  v ,  >  >  0  >  ti(v:) 
and  I  /(  u (  v  , ) )  / ( ti (  v , ) )  (  g  Dr'  "  for  some  positive  constant  D. 

Proof:  Since  u\x)  >  0  for  .v„  <  a  g  1.  and  <p(.v„)  =  -  c( .v,,)/i(  u{  v(, ))  >  0 
and  a  "(a,,)  =■  f  '<j  ( .v„).  it  follows  from  (i)  of  Lemma  3.6  that  there  exists 
A.,  a,  >  ,v,,.  with  <r(  v  i )  =  0,f'  4  and  |w(.v(l)  -  k(a:)|  +  |.v,  -  x0|  =  0(1)f  s. 
Moreover,  u (  v , )  <  -O  for  some  positive  0  independent  of  f.  Similarly,  front 
(ii)  of  Lemma  3.6.  there  exists  a,  with  a,  <  a„  and  (  v , )  =  -OiF1  4.  w ( a , )  >  C 
and  1  a,  -  a ,,|  =  0(1  )«■■'  v.  (t  remains  to  verify  the  last  estimate  in  the  lemma. 
Integrate  (1.1)  from  a,  to  v(l  to  obtain 

Ft/'(.\,)  -  ft/'Lx,)  =  /(u(.\,))  -  /  (  M  (  A ,  )  )  -  j  c(  x)h(ll(  \))  d\ . 


Since  clt  =  0(1).  ( v,  -  v,|  =  0(1)f,/I<  and 


u'(  a  )  = 


<P(A)  +  c(  A  )  ft  (  l<  (  A  )  ) 

/>(*)) 


which  is  hounded  at  v  =  a,,  a,  because  |u(.v)|>  O  there,  we  have 

/(«(*:))  -  /(«(  v,)>  =  0(1  )f  +  0(  1 ) ( a 2  -  a,)  =  0(\)r'  \ 


This  proves  the  lemma. 
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Theorem  3.8.  Lei  u(x)  =  m(x.  f)  be  a  solution  of  (1.1),  (1.2).  Then  for  t 
sufficiently  small.  w(x,  f )  is  close  to  an  inviscid  time-asymptotic  solution  of  (2.1) 
with  the  same  boundary  data  u ,  and  ur.  More  precisely,  there  exists  an  inviscid 
solution  u(  x.  0)  with  boundary  data  uf  and  ur  such  that 

(i)  if  u(  x,  0)  has  a  boundary  layer  («,.«,)  at  x  -  0.  then  there  exists 
0  <  x,  =  0(l)f  ’/f<  such  that  u(x, )  -  tq  =  0(  l)f3/R  and  u(x)  is  monotone  oner 
(0..v, ); 

(ii)  if  u(x.  0)  has  a  standing  shock  wave  (  u2.  u})  at  x  =  x0.  then  there  exist 
x,.  x3.0  <  x,  <  .v,  <  1.  such  that  |«(x2)  —  u2|  +  |m(x3)  —  u3|  =  0(l)e3/R, 
|x3  —  x,|  =  0(l)e3  R  and  u'(x)  <  0  for  x2  <  x  <  xy, 

(iii)  if  u( x. 0)  has  a  boundary  layer  (u4,  ur)  at  x=  1.  then  there  exists 
x4. 1  —  0(l)r/R  <  x4  <  1.  such  that  i/(x4)  —  «4  =  0(1)£3,,R  and  u(x)  is  mono¬ 
tone  over  (,v4.  1 ). 

(iv)  outside  all  possible  layers  of  the  types  (i)— (iii). 

|/(m(x))'  -  c(x)h(u(x))\  =  0(  1  )£1/4  and  |m(x)  -  u(x.O)!  =  0(1  )r  *. 


Proof:  From  Proposition  3.5  we  see  that  a  monotone  solution  (u(.x  )  of  (1.1). 
(1.2)  is  close  to  being  inviscid  except  for  possible  boundary  layers.  Lemma  3.6 
says  that  the  boundary  layers  almost  satisfy  the  inviscid  boundary  conditions 
(2.4)  and  (2.5).  Consequently,  we  can  find  u,.  5,.  fi4  and  ur  with  the  property  that 
( u,.  m, )  is  an  admissible  inviscid  boundary  layer  at  x  =  0.  (u,.  u4)  an  inviscid 
stationary  wave.  (u4.  ur)  an  admissible  inviscid  boundary  layer  at  x=  1,  and 
I u,  -  M/l  +  |  Si  -  m(xi)  |  +  |u4  -  «(x4)  (  +  \ur  -  ur  1  =  0(1  )ev*.  for  some  x,.  x4 
with  the  prescribed  properties  in  (i).  (iii).  (iv)  of  the  theorem.  We  denote  by  v(x) 
the  inviscid  time-asymptotic  state  with  boundary  values  u,  and  ur.  From  the 
inviscid  theory  of  Section  2,  there  are  at  most  three  solutions  to  the  inviscid 
problem  with  given  end  states.  Moreover,  the  inviscid  stationary  solutions  and 
the  layers  depend  smoothly  on  its  boundary  values.  Since  | u,  -  u,\  +  \ur  -  ur |  = 
0(l)e3/x,  we  conclude  that  there  exists  an  inviscid  solution  (w(x,0)  with  bound¬ 
ary  values  n,  and  ur  such  that  u(x, 0)  is  close  to  u(x)  in  the  sense  of  (i)-(iv). 

For  a  nonmonotone  solution  w(x),  it  follows  from  Lemma  3.1  that  u(x)  has  a 
unique  minimum  at  x  =  x0.  There  are  two  cases:  (I)  x0  is  close  to  x  =  0  or 
x  =  1;  (II)  |x0|  >  Cf  and  |x0  -  1|  >  Ce3/8  for  some  large  positive  constant 
C.  In  the  case  where  x0  is  close  to  x  =  0.  we  have  «'(x)  <  0  for  0  s  x  g  x0. 
u'(x0)  =  0  and  so.  by  integrating  (1.1)  from  0  to  x3. 


ru'(x(l)  -  eu'(O)  =/(«(x0))  -  f(u(0 ))  +  f  "c(x  )h(u(x))  dx, 

*'() 


-/(“(*o))  +  /(“( 0))  <  (  °c(x)h(u(x))  dx  =  0(1  )x0  =  0(l)eVR. 

•'O 
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Thus  u(x)  has  a  boundary  layer  near  x  =  0  which  satisfies  the  inviscid  boundary 
condition  (2.4)  except  for  a  small  error  0(l)c3/8.  Similarly,  if  x0  is  close  to  1, 
then  u(x )  has  an  admissible  boundary  layer  at  x  =  1.  Thus  the  above  arguments 
for  monotone  solutions  can  be  applied  to  u(x)  outside  the  layers  and  Case  (I)  is 
treated  accordingly.  For  Case  (II)  we  apply  Lemma  3.7  to  locate  x2  and  x3  close 
to  x0,  x,  <  x0  <  x3,  so  that  u(x)  has  an  interior  layer  in  (x2,  x3).  The  solution 
h(x)  is  monotone  over  (0,  x2)  and  over  (x3, 1)  and  so  the  above  arguments  for 
monotone  solutions  apply  again.  This  completes  the  proof  of  the  theorem. 

Remark  3.9.  The  above  theorem  does  not  yield  an  optimal  thickness  of  the 
layers  and  the  distance  between  viscous  and  inviscid  solutions  outside  the  layers. 
Suppose,  for  instance,  that  u(x)  has  a  boundary  layer  at  x  =  0  and  that  ut  >  0. 
From  Theorem  3.8,  it  is  a  shock  layer,  u'(x)  <0  for  0  5  r  <  x0,  u'(*o)  = 
0,  u(x0)  <  0  for  some  x0  >  0.  Moreover,  for  e  sufficiently  small,  f(u,)  < 
/(u(x0)).  From  (1.1)  we  have 

u"(x)  =  7(/'(“U))u'(*)  -  c(x)/t(u(x))). 


eu'(x)  =  f(u(x))  -  f(u(x 0))  +  (  °c(y)h(u(y))  dy,  0  g  x  g  x0. 

J  X 

From  the  first  identity  and  (1.5),  ch  <  0,  we  see  that  u"(x)  >  0  for  x  g  x  g 
x0,  u(x)  =  0.  Choose  x,.  x  <  x3  <  x0,  with 


f(u(xn))  -/(u(x,))  =  j  °c(y)h{u{y))  dy, 
J  X\ 


or 

|u(x0)  -  U  (  X|  )  j  =  0(l)|xo  -  X  j  | , 
for  some  positive  bounded  0(1),  Since  «'(x0)  =  0,  we  have 

\f(u(x))  -  f(u(x0))\ <  f  °c(y)h(u(y))  dy  for  x,  <  x  <  x0, 

J X 

and  so 

eu'(x)  =  0(l)(x  -  x0), 

for  some  bounded  and  positive  0(1)  and  x,  <  x  <  x0.  Integrate  to  yield 


|m(x0)  -  «(x,)j  =  ^7^-|x0  -  x,|2. 


£ 
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Thus  we  have 


0(1) 


l-*o  -  *il2  =  O(l)|x0  -  x,|. 


or 


1*0  -  *il  =  0(l)f,  !«(x0)  -  u(x,)|  =  0(l)e. 

Between  x  and  x,.  we  have 


|/(u(x))  -  /(u(x0))|  >  f  °c(y)h(u(y))  dy , 

J  X 


SO 


£li'(x)  =  0(l)(/(u(x))  ~/(u(x0)),  X  <  X  <  X,. 


Integrate  it  from  x  to  x^ 


/•o  _ du 

V<,i/M*o))  -/(«) 


=  0(l)(x,  -  x). 


or 


0(1  )e log e  =  0(l)(x,  -  x). 

Thus  we  conclude  that  x0  -  x  =  0(l)eloge,  a  similar  estimate  holding  for 
x  =  0.  Thus  the  thickness  of  the  boundary  layer  is  of  the  order  e  log  e.  The  same 
holds  for  other  layers  as  well.  Details  are  omitted. 

4.  Local  Uniqueness  and  Bifurcation 

We  want  to  establish  a  local  uniqueness  theorem  which,  when  combined  with 
the  a  priori  estimates  of  the  last  section,  determines  the  exact  number  of  solutions 
of  (1.1),  (1.2).  For  this  we  employ  the  shooting  method.  Let  u(x)  =  u(x ,  e,  jB)  be 
a  solution  of  (1.1),  with  initial  slope  /3: 

(4.1)  u(0)  =  u,,  u'(  0)  =  0. 

The  following  crucial  lemmas  establish  the  dependence  of  u(  x ),  x  >  0,  on  the 
initial  slope  /?.  Set 


w(x)  =1 


<M*) 

dp 


(4.2) 
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It  follows  from  differentiating  (1.1)  and  using  (4.1)  that 
ew"  —  f'(u)w'  +  mw  =  0. 

(J  m  =  m(x)  =  c(x)h'(u(x))  -  f"(u(x))u'(x). 

(4.4)  w(0)  =  0.  h'(0)  =  1. 

It  is  often  convenient  to  rewrite  (4.3)  as 

(4.5)  fit'"  -  ( f(u)w )’  +  c(x)h'(u) h-  =  0. 

The  first  lemma  deals  with  the  easy  stability  case.  (1.6). 

Limma  4.1.  Suppo.se  that  e(x)h’(u)  <  0.  Then  w(x)  >  0  for  all  0  <  x  <  1. 
Proof:  Integrate  (4.5)  repeatedly  to  yield,  for  U  ,v  >  x'  g  0. 

Fie  (.x)  =  rif'(.v')  +  /'(«(: v  ))»•(. tr)  -  f'(u(x"))  ii(.v')  -  J  ch’w(r)  dr. 


,  ,  t  '  \  I  f'f'(u)U)  ,r 
it  ( -V  )  =  it  (  V  )exp  |  - - - di 


M  /'(  m(  -v'))  /  M 

H  (  A  )  —  — - - -  If  {  A  ) 


'  i  rnytt)) 


Hr  l 


dt  dv 


7  /  dyf  e(T)h'(  u(  t))w(t  )exp  ^  -  j  — —-d^dr. 

Setting  v'  =  0.  it  follows  since,  by  hypothesis,  ch'  <  0  that  iv(.v)  >  0  so  long  as 
m  (  i  )  >  0  for  0  <  y  <  v.  By  (4.4)  it  is  clear  that  if (  v )  >  0  for  y  close  to  0.  This 
proves  the  lemma. 


Theorem  4.2.  Suppose  that  c(x)h'{u)  <  0.  Then  (1.1).  (1.2)  have  a  unique 
solution  w  hich  tends  to  the  corresponding  inviscid  solution  as  e  — ►  0  , . 

Proof:  Uniqueness  follows  from  Lemma  4.1  or  by  the  maximum  principle. 
The  existence  of  solutions  is  established  by  the  shooting  method.  We  defer  this 
until  later  when  we  deal  with  the  issue  for  the  unstable  case  in  Theorem  4.7. 


For  the  remainder  of  this  section  we  treat  the  more  interesting  and  much 
harder  instability  case  e(x)hfu)  >  0.  For  this  we  need  to  look  at  the  boundary 
layers  and  interior  layers  separately.  The  first  lemma  on  a  single  boundary  layer 
refines  the  classical  result  of  Coddington  and  Levinson  (2).  Throughout,  we 
assume  that  r  is  small. 
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Lemma  4.3.  Suppose  that  u(x)  is  a  solution  of  (1.1).  (1.2)  and  has  only  a 
boundary  layer  at  x  —  0.  That  is.  there  exists  xt)  =  0(1  )e  such  that  (2.4)  holds  with 
w()  =  u(x0)  andu'(x)  —  0(1)  independent  of  e  for  x  >  x„.  Then  the  solution  h(x) 
of  (4.3).  (4.4)  satisfies  w(x)  >  O  for  1  2  .?  >  0  and  w''(x)  >  —  \e  'uLv)  for 
x  >  Del/1  and  a  constant  \  —  Cel/2,  where  C  and  D  are  positive  constants 
independent  of  e. 

Proof:  We  first  show  that,  for  some  x,  =  0(1)e1/:. 


(4.6) 


w(x{)  >  Ce,  m'(.v,)  =  0(1  )f1/J, 
u'(x)  =  0(1),  jc,  g  x  g  1. 


The  latter  estimate  follows  from  the  results  of  the  last  section.  Remark  3.9. 
since  v,  lies  outside  the  layer.  Our  proof  for  (4.6)  is  carried  out  along  the  same 
line  as  similar  arguments  in  Coddington-Levinson  [2j.  except  for  some  refinement 
due  to  the  boundary  layer.  We  note  that  in  [2]  the  authors  dealt  with  only  the 
subsonic  case  and  we  need  to  treat  the  transonic  case.  Integrate  (4.5)  from  0  to  v 
to  yield 


(4.7),  fit ■'(  a  )  —  /'(  u  )  ix(  x  )  =  f  -  f  (  (  r  )/('( t/(  r  ) )  u(  r  )  dr. 

•'ll 

Another  integration  results  in  an  expression  for  u(.x ): 

(4.7),  u  (  v)  -  /  /•.  (  v.  s  )  f/v  -■  -  (  c{r)h'(u(r))w(r)  I  li{\.s)dsdr, 

•'ll  *  -Yi  J , 

where 

/-.  (  x .  V  )  =  exp  j  '  j  /'(  u(  t  ))  dr  |  . 

Let  £  =  0(1  If1  :  he  such  that  f'(u(  a  ) )  g  -  A  <  0  for  £  s  x  <  1.  for  some 
A  >  0  and  V  =  max' |/'(«(  a  ))|:  0  <  v  g  £).  Set  -  £  +  l.f  \  where  I.  - 
1  +■  ()( I  )  .VA  '.  Then,  for  Os.vg  £. 

IA  £.  x)  <  exp  j  -  (4  v)  ■  g  exp{  MO{  1  )r  1  :  j . 

TA  x.  A)  1  exp  J  ~  (  v  --  £ )[  g  exp{  -  kl.t  1  : } . 

/(v.  v)  <.  exp{  (  -  A7.  +  NO(  1))f‘  : }  §  exp(  -Af  1  :}. 
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Thus,  for  x  g  X|, 

Ce(x,s)cIs=  f(E(x,s)ds+  fXE{x,s)ds 
J0  Jo  £ 

g  J£exp{ ds  +  jTexpj  -  j(x  -  s)}  ds 

2  £ 

^  T- 

Similarly,  the  second  term  in  the  right-hand  side  of  (4.7)2  can  be  written  as  the 
sum  of  /,  and  I2,  where 

;  IfV(r)fi'(u(r))w(r)[/££(x,s)^+  f E(x,s)ds  dr , 

1  £  Jo  r  £ 

1 2  =  -  j^'Vc(r)A'(«(r))w(r)jT,£(x.5)d«/r. 

If, |  S  J^|w(s)![exp{-fce'1/2i}  +  I ]  =  X  ^M*)!** 

|/2|  g  ~  jfV(s)|<fe,  where  k(x)fi'(«(x))|  g  n.  Ogxgl 

Combining  these  estimates  we  obtain 

|w(x)|  Sy  +  j|  M*)l* 

and  this  in  turn  implies 

k(x)|  S  Texp{¥} 

and 


w(x)  «  fXE(,x,s)ds  +  xO(t). 


Thus 
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As  in  [2],  we  may  show  that  by  taking  care  of  the  boundary  layer, 


w(x)  = 


[l  -  £(*.£)]  -  77 


o(t 


xO(e) 


-/'(“(*)) 
for  x  s  x,,  and  by  (4.7),  we  have 

jw'(x)  j  g  1  +  O(e)  +  X0(1)  + 


f'(u(x))  f'(u(x)) 


exp 


2  nx 


-  1 


or 


w'(x)  =  0{x  +  e)  for  x  >  x,. 

This  proves  (4.6). 

Let  X  be  a  small  positive  constant  to  be  determined  later  and  set 

Xx 


We  have 


w(x)  =  n(x  )expj  j-  j 
w'(x)  =  ju(x)  +  t/(x)expj  - 


From  (4.6)  we  have 


c(x, )  >  0,  c'(X|)  >  0. 


provided  that  X  s  Ce1  11  for  some  large  C.  By  (4.3)  and  the  above  it  follows  that 
ev"  +  (  -/'( u )  +  2X  )  c' 


(4.8) 


X2  -Xf'(u) 


+  f"(u)u‘  +  c(x)h'(u) 


1-  =  0. 


In  view  of  (4.7)  we  know  that 

\f"(u)u'  +  o(x)fT(«)|  =  0(1). 

for  x,  g  x  g  1.  By  Theorem  3.8,  outside  the  layer,  x  >  x,.  u(x)  is  close  to  the 
invisctd  stationary  solution  which  is  subsonic,  f'(u(x))  <  0.  Thus  we  can  choose 
X  small  such  that  the  bracket  in  (4.8)  is  positive.  One  such  choice  is  X  =  Ce1  2  for 
some  large  C.  Thus  the  lemma  is  proved  by  the  maximum  principle, 

Lemma  4.4.  Suppose  that  u(x)  is  a  solution  of  (1.1),  (1.2)  and  has  only  one 
boundary  layer  at  x  =  0  or  x  -  1,  but  not  both,  and  has  no  interior  layer.  Then  the 
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solution  »  (  v  )  of  ( 4 . 3 >  with  either  »•(())  =  0.  vc'(0)  =  1  or  u(l)  =  0.  u  '( 1 )  =  1 
sattstie i  either  » ( I )  >  0.  » '<  1  (  >  -\f  lu(l)  or  >i(0)  <  0.  vx'(0)  <  \i  'n(O)  for 
some  A  =  (  i 1  . 

Proof  Suppose  that  u  (0)  =  0.  vx  ’(0)  =  1.  The  lemma  is  reduced  to  the 
previous  one  if  the  boundary  layer  is  located  at  x  =  0.  Suppose  that  the 
houndarx  laver  is  located  at  v  =  1.  Then  by  the  same  argument  as  in  the  last 
proof  xv e  can  show  that  the  solution  if(.v)  of  (4.3)  with  tc(l)  =  0.  vc'(l)  =  1 
satisfies  iT ( ( ) >  <  0.  Since  the  Wronskian  vc(  v)iv'(  x)  -  iff x )»••’(. v)  does  not  change 
sign  and  is  positive  at  \  -  o.  we  see  that  u  <  I )  >  0.  Other  parts  of  the  lemma  are 
proved  similar!). 

l.i  MM  a  4.5.  Sttpi  >ose  that  tt(  x  )  is  a  solution  of  ( 1 .1 ).  ( 1 .2)  with  only  an  interior 
laver.  that  is.  there  exists  x„  and  x". 0  <  ,v„  <  a1’  <  1.  such  that  j//'(.v)|=  0(1) 
independent  of  t  for  .x  £  (  v,,.  A1').  |  v1' -  x0|  =  0(  1  )e. /(  «( A11))  =  ./(  u(  v,, ))  + 
Ol  1  )t.  at  \ ,, )  >0  -  ii(  v").  77ie/i  the  solution  h(a)  of  (4.3).  (4.4)  satisfies  ic(.v)  >0 
for  \  <  x  ,.  u  (  \  j )  r-  I).  » (  x )  >  0  for  x  >  .v,  for  some  .v,  not  in  the  interior  laver. 
x1  -  x , ,  with  it,  -  tnj  -  O())r|)ogr|.  Moreover.  h'(Aj)  <  0.  and/or  some  con¬ 
stants  (  .  I)  ().»>■'( .v)  <  Af  'ti(.v)  or  x  >  .v,  +  /V  1  and  \  =  tV  :. 

Proof:  The  proof  consists  of  three  steps  investigating  the  behavior  of  n  (.v) 
before  the  layer,  in  the  layer,  and  after  the  layer. 

Step  I.  l  et  A  be  a  positive  constant  to  be  determined  later  and  set 
(4.4)  h  (  v  )  =  r(  v  )exp|  7-  j  . 


We  have  from  (4.3).  (44)  that 
14.10) 


w'{  x )  -■*  *  »  (  v  )  +  e'(  x  )exp 


(4.1 1  ) 


h"  a  (2A  -  f(u))r‘  + 

r(0)  =  0.  e'(0)  =  1. 


A-  ,  A 

7  -  /  ( w )  7  +  m 


r  -  0. 


Prom  t  heorem  3.H  we  know  that  u(x )  is  close  to  a  supersonic  inviscid  stationary 
wave  for  0  <  x  <;  \ Thus  we  may  choose  a  small  A  such  that  (A:  - 
ft  m(  v  ))A  )r  1  «  -  1  for  0  <  x  g  a0.  Since  mix)  is  bounded  over  0  g  v  §  ,v„. 


T  “/'(«)  7  +  m  ^  ~ 


0  <  A  <  A,,. 


Thus,  by  the  maximum  principle,  applied  to  (4.11).  (4.9)-(4.11)  imply  that 
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r(  a  )  >  0.  r’(  v )  >  0  and 


«■'  >  —  if  for  0  <  x  <  x„. 


A  f'(u) 


This  and  (4.3)  yield 

f«"  =  /'(  u )  if'  —  m\ f  > 
and  so.  for  small  t. 

Fif "  >  w'  for  0  S  r  |  x„. 

This  and  ihe  iniiial  data  (4  4)  for  it  yield 


w  +  jf'(u)w'. 


if '(  a  )  >  exp 


(4.12) 


I  *i-M 

\  r 


"lx)>  ^(exp!^r)  - 1)- 


0  S  A  S  X„. 


for  some  A,  >  0. 

Step  2.  Define  a*  by  u(x*)  =  0.  v0  <  x*  <  a".  Integrate  (4.5)  front  0  to  a*. 
to  obtain 


f*e'(.V*)  = 


=  f  ~  ('  f(  x  )h'(  u(  a))  if  (a  )  dx. 

•'ll 


Since  x*  is  within  the  layer  we  have  from  Lemma  4.4  that  ie(.x)  >  0  for 
0  <  a  s  a*.  Moreover,  we  have  the  estimate  (4.12).  Thus  the  above  yields 
it '(a*)  <  0  and  we  conclude  that 


(4.13) 


w '(a)  =  0  for  some  x,  x0  <  x  <  x*. 

k  '(x)  <  0  for  x  <  x  s  x*. 


Step  3.  We  want  to  show  that  there  exists  x,  >  x.  |xi  -  x|  small,  such  that 
h-(x,  )  =  0.  For  the  moment  we  assume  that 


(4.14) 


w(x)  >0  for  0  <  x  <  x,. 


for  some  x,  >  x.  We  have  either  x,  =  1  or  x,  <  1  and  w(x,)  =  0.  From  Lemma 
4.4  we  know  that  x,  does  not  lie  in  the  interior  layer.  Integrate  (4.5)  from  x  to 


54 


S.-B.  HSU  AND  T.-P.  LIU 


.x,  x  <  x  <  x ,,  then 


(4.15) 


ew'(x)  =  f'(u(x))w(x)  -  f'(u(x))w(x) 

~  (  c{y)h’{u(y))w(y)  dy. 

Ji 


From  (4.14),  (4.15)  and  u(x)  >  0  we  obtain 

(4.16)  ew'(x)  <  f'(u(x))w(x)  for  x  <  x  <  xL. 


Then  (4.13).  (4.16)  yield 

(4.17)  w'(x)  <  0  for  x  <  x  <  xt. 

Let  ,x  =  x  +  Del/4  for  some  large  positive  constant  D.  Then  x  is  outside  the 
layer  and  x  >  x.  where  .x  is  characterized  by  u'(x)  ~  0.  Suppose  (4.14)  holds  for 
x,  >  .x  +  Z5 ( e 1 1  og  f  |  for  some  large  positive  constant  D j.  Integrate  (4.3)  from  x 
to  x: 


(4.18) 


w'(  .x ) £ ( .x )  =  w'{ x )  —  —  (  m(s) E(s)w(s)  ds , 

e  Jx 

£(.x)  =  exp |  -  7 /_/'(«(>’))  4'j- 
m(x)  =  c(x)h’(u(x))  -  f"(u(x))u'(x). 


Divide  (4.18)  by  £(.x)  and  integrate  to  obtain 


w(.x)  =  w(x)  +  w'(x)J  E(y)  1  dy 


(4.19) 


~jj_E(y)  1  dyj_m{s)E(s)w(s)  ds. 


From  (4.14),  (4.15),  (4.17)  we  obtain 


1 


*>'(*)  <  7 (/'("(*))  -  0(l)(x  -  x))w(x). 


Thus  (4.19)  and  the  above  yield 


w(x)  <  w(x) 


(4.20) 


i  +  fW)) n  j;£(  v)-> 


-  7  f_  E(y)~ldy  j^m{s)E(s)w(s)  ds , 
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for  x  ^  x  g  Xy.  Integrate  (4.16)  to  obtain 


j(s)  5  w(x)expj^ 


di\  for  rS  j$  x^ 


Since  m(s)  is  bounded,  the  above  estimate  yields 

|m(s)£(s)w(s)|  ^  Cw(i)  for  x  ^  x  <;  xt 
for  some  C  independent  of  t.  This  and  (4.20)  imply 

/'(“(*))  -  £>(l)(x  -  x)  r' 


H'(x  )  <  H'(x) 


(4.21) 


1  + 


I  E  (  v )  ‘  dy 


+  E(y)  \y  -  x)  dy 

Take  x  =  x  +  2e  |  log  e|  in  (4.21).  Then  we  obtain 


w  ( x )  <  w  ( x ) 


1  +  "  U£  X~  /-  E(y)  1  dy  +  E{  v)' '  dy 


(4.22) 


X(0(l)2f|loge|  -  0(1  )e1/4 


<  w(x) 


1  + 


/'("(x)) 


/  E(y)  ~'  dy 


for  e  >  0  sufficiently  small.  It  remains  to  show  that  the  above  bracket  is  negative 
and  thus  contradicts  (4.14).  For  this  we  expand  E(y)~ 1  by  Taylor  expansion, 
taking  note  of  the  fact  that  u’(x)  =  0,  x  is  in  the  interior  layer,  and  x  >  x  is 
outside  the  layer, 

(4.23)  /'(«(*))  =/'(“(*))  +  (*- •*)«'(£)/"(“(?)).  x  <  f  <  x. 

Since  |  <p  ( x )  |  g  Clei/A  for  x  2  x°, 

_  ~C(x)h(u{x))  -  <p(x) 

U(X)~  -/'(u(x))  -  2’ 


for  some  positive  C2  and  x  >  x.  Thus  we  have  from  (4.23)  and  (1.3),  f"  >  0,  that 
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for  some  C  independent  of  e. 

Consequently, 

FA  i 

)  1  =  cxp|  7  /  /'("(  ' 

OMxj 

(4.24) 

"^!ij 

v))  +  C(x  -  x))dxj 

=  vxp[  -  ~\  r(u( 

a  ))  |(  r  -  x )  j-expj  ^(y  -  x)Jj.  x  <  y  <  a,. 

Let  SI 

■■  /’(  u(  V  ))  . 

Then 

uf'i: 

(  r  )  '  dv  >  SI 

l'exp\ 

-~(y  ~  a  )|e.\pj  yi  <  >'  -  *r) 

>  si 

l'c*  ?! 

■-  ~(.V ■  -  -v)j(l  +•  Yt  <  .1  -  V )-)  dy 

—  t 

I  exp 

I  SI  .  _,i  Cr!  r\U\  si,  , 

-  (  a  x )  +  -  -  ,  /  :  e  d:  >  e 

1  f  1  j  2  SI  -  Jo 

Thu'  the  bracket  in  (4.22)  is  negative. 

To  finish  the  proof  we  observe  from  (4.15)  and  w(.v,)  =  0  that  u'(x,)  ^  0. 
Since  v,  is  on  the  right  side  of  the  interior  layer,  we  may  apply  Lemma  4.3  to 
-  vv(  v  (/«’(  v,)  to  show  that  >r(.v)  <  0  for  A  >  v,. 

l.i  mm  a  4.6.  Suppose  that  u(x)  is  a  solution  of  (1.1)  which  either  luis  tin 
interior  lover  not  located  within  Of]  )t  of  x  =  0  or  x  =  1,  or  has  no  interior  lover. 
In  the  former  ease,  the  corresponding  solution  u(  a)  oj  (4.3).  (4.4)  has  the  property 
that  u  ( ! )  <:  0.  and  m  the  latter  ease,  w  ( 1 )  >  0. 

Proof:  Using  l  emmas  4.3  4.5  and  Theorem  3.8,  it  remains  to  treat  the  ease 
where  u(x)  has  either  an  interior  layer  and  a  rarefaction  layer  at  v  =  1.  ef.  (ii)  of 
Case  HI  in  Section  2.  or  a  houndary  layer  at  x  =  0  and  a  rarefaction  layer  at 
v  -  1.  cf.  (ni)  of  Case  F.l  and  Case  D1  in  Section  2.  In  all  these  cases  the  solution 
is  subsonic  before  the  rarefaction  layer,  which  connects  the  sonic  slate  zero  to  the 
positive  stale  u,  around  x  =  1.  We  treat  the  case  without  interior  layer;  the  one 
with  the  interior  laver  is  similar.  From  Lemma  4.3  we  have,  before  the  rarefaction 
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Since  C,  and  8  are  positive  constants  independent  of  e  and  X  =  CeI/2,  the  above 
estimate  yields,  for  small  e, 

f  <Sj(-m(s))£(i)w(i)  ds  »  — 

x0 

This.  (4.25)  and  (4.28)  imply 

h'(1  —  5)£(1  —  8)  >  — Ae_1w(x0)  +  J  —  (  —  m(s))E(s) w(s)  ds  >  0. 

x0 

Thus  we  have 

(4.29)  w'(x)  >  0  for  x  =  \  -  8. 

This  and  (4.25),  w(x)  >  0  for  x  =  1  -  S,  implies  that  w'(x)  >  0  for  1  —  8  g  x 
<  1.  Indeed,  from  (4.3)  if  w'fx)  =  0  for  some  x.  1  —  6  <  x  g  1  and  w(x)  >  0. 
then 


fu"(x)  =  -m(x)w(x), 

which  is  positive  by  (4.26),  a  contradiction.  Thus  w'(x)  >  0  for  1  -  8  <  x  g  1. 
In  particular.  w(l)  >  u  (l  -  8)  >  0.  This  proves  the  lemma. 

Theorem  4.7.  Suppose  that  c(x)h'(u)  >  0.  Then,  for  sufficiently  small  t. 
(1.1).  (1.2)  have  one  or  three  solutions  which  tend  to  the  corresponding  inriscid 
solutions  as  e  -»  0,.  Moreover ,  (1.1),  (1.2)  have  three  solutions  for  small  e  in  Case 
El  and  (iii)  of  Case  E3  (in  Section  2  for  inviscid  classification),  and  in  other  cases 
there  exists  only  one  solution. 

Proof:  That  there  exist  at  most  one  or  three  solutions  in  each  respective  case 
follows  from  the  a  priori  estimate  in  Section  3  in  linking  the  solutions  to  inviscid 
solutions  and  the  local  uniqueness  theorem  as  a  consequence  of  Lemmas  4.4-4.6. 
It  remains  to  show  the  existence  of  a  solution  which  is  close  to  any  given  inviscid 
solution.  This  is  done  by  the  shooting  method.  We  shall  carry  out  the  analysis  for 
the  case  where  there  is  a  shock  layer  at  x  =  0  and  a  rarefaction  layer  at  v  =  I : 
other  cases  can  be  treated  by  similar  arguments.  For  simplicity,  we  assume  that 
«,  >  u,  >  0  and  ur  >  0.  The  inviscid  solution  consists  of  a  shock  layer  ( uh  «,)  at 
x  =  0,  an  inviscid  stationary  solution  (u,.0)  satisfying  /(u),  =  c(x)h(u)  and  a 
rarefaction  layer  (0.  ur)  at  x  =  1.  We  want  to  find  a  value  a  =  «0  so  that  the 
solution  u(x.  a)  of  (1.1).  with 

u(0,a)  =  u  /  •  u'(0,a)  =  a. 

satisfies  u(l.a(l)  =  ur.  For  |a|  »  l,a  <  0.  i.e.,  /’(u,)a  -  c(0)h(u,)  g  —  C,f*  4. 
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we  have  from  (ii)  of  Lemma  3.4,  that  u(x,  a)  has  a  boundary  layer  at  x  =  0. 
From  (1.1),  u"  >  0  so  long  as  u  <  0  and  u'  <  0.  Thus  u'(x0,  a)  =  0  for  some 
x0  >  0.  By  Remark  3.9  on  the  thickness  of  boundary  layer,  we  have  x0  =  0(l)e. 
Integrate  (1.1)  over  0  g  x  g  x0  to  get 

-ea  =/(u(x0,a))  -  f(u,)  -  (  °c(x)h(u(x,a))  dx , 

Jo 


or 


f(u(x0,a))  +  0(l)emax{/i(u)  :  u(x0,a)  <  u  <  u,}  =  e|a|  +/(«,). 

Thus  from  (2.7)  we  see  that,  as  a  -»  —  ao,  u(x0,a)  -*  —  oo.  For  the  boundary 
layer  to  exist  at  x  =  0  it  is  sufficient  to  have  /'(«/)«  ^  c(0)/i(u,)  -  C1e1/4.  For 
a  =  a,  satisfying  f'{ul)al  =  c(0)h(u,)  -  C)£1/4  the  above  identity  yields 

/(«(x0.a,))  =/(«,)  -  eai  +  0(l)e 

-/(«/)  -  e(r(u,)y'{c(0)h(u,)  -  C,e1/4)  +  0(1)« 

=  /(«/)  +  0(l)e. 

Thus,  as  f  -*  0,  u(x0,  n,)  -»  u,>  and  so  for  small  e,  u(x0.  a,)  >  u,.  Since 
u(x,  a)  is  close  to  being  inviscid  for  x  >  x0  so  long  as  u(x,  a)  is  not  close  to  the 
sonic  state  zero,  (i)  of  Lemma  3.4.  it  follows  that,  for  small  c,  u(x.a)  stays 
subsonic  for  u(x0,  a)  <  u,  and  becomes  sonic  at  x  =  x(a)  <  1  for  u(x0,  a)  > 
u*.  The  latter  holds  for  a  =  a,.  Moreover,  once  u(x(a),  a)  =  0  for  x(a)  <  1, 
which  is  the  case  for  a  =  a,,  u(x,  a)  has  a  boundary  layer  at  x  =  1,  Theorem  3.3. 
and  is  strictly  increasing  for  x0  <  x  <  1.  Moreover,  from  Lemma  3.5,  if  1  -  x(a) 
>  0( l)el/\  then  u(x,  a)  becomes  +  oo  before  reaching  x  =  1.  Since  u(l,  a)  is  a 
strictly  increasing  function  of  a  for  a  <  a,.  Lemma  4.3,  and  0(1.0,)  = 
x.  w(l.  a)  -*  -  *  as  a  -*  -  oo,  we  see  that  there  exists  an  a0  <  a,  with  u(  1,  a0) 
=  ur.  Moreover,  since  a0  <  at,  u(x,  a0)  has  a  shock  layer  at  x  =  0.  By  Theorem 
3.8.  «(x,  a0)  also  has  a  rarefaction  layer  x  =  1.  We  have  thus  constructed  the 
designated  solution  of  (1.1).  The  proof  of  the  theorem  is  complete. 

5.  Asymptotic  Stability  and  Instability 

In  this  section,  we  study  nonlinear  stability  and  instability  of  stationary 
solutions  of 


+/("),  =  cuxx  +  r(x)Mu). 
u(0.  t\=u,,  u(I,f)  =  u,.  m(x,0)  =  mh(x). 


0  g  x  g  1 . 


(5.1) 
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Let  U{x)  =  U{x,  e )  be  a  stationary'  solution. 


(5.2) 


f(l')x  =  eU„  +  c(x)h(U), 
V(0)  =  u,.  {/(!)  =  ur. 


Set  n(  v,  t)  =  u(x.  t )  -  U(x).  From  (5.1).  (5.2)  we  have 

w,  +  /'( w  +  U  )(uv  +  U\ )  =  fH-„  +  f’(U)(Jx-  c(x)h(U)  +  c(x)h(  tr  +  Li). 
(5.3)  v.(0, /)  =  0,  vi'(l.r)  —  0. 


The  linearized  equation  is 

»c,  +  /'(  U  )  wx  +  f"[U)Uxw  =  £H\(  +  c(x)h'(U  )  w. 

H'(0.  t)=0.  H'(l./)=0. 

We  shall  use  spectral  analysis.  Set  u(.v,  /)  =  ex'q(x)  and  obtain  from  (5.4) 


tq"  -  (f'(U)q)'  +  c(x)h'(U)q  =  At). 

4(0)  =0.  4(1)  =0. 


Thkorlm  5.1.  Suppose  (1.6)  holds.  c(x)h'(u)  <  0.  then  every  steady  state 
L(x)  of  (5.1)  is  asymptotically  stable.  When  (1.7)  holds.  c(x)h’(u)  >  0.  then 
steady  states  containing  no  interior  layer  are  asymptotically  stable:  those  containing 
interior  layer  are  asymptotically  unstable. 

Proof:  A  steady  state  U(  x )  is  stable  if  the  largest  eigenvalue  A  is  negative. 
From  linear  Sturni-Liouville  theory  (see  [1))  the  eigenfunction  q(x)  correspond¬ 
ing  to  the  largest  eigenvalue  A  is  of  one  sign.  We  may  therefore  assume  that 


(5.6)  4 (  v  )  >  0.  0  <  ,x  <  1 .  4'(0)  =  1 .  4'(1)<0. 

When  <-(.v )/»'(«)  <  0.  we  have,  integrating  (5.5).  that 

(5.7)  f(4'(l)  -  4'(0))  +  f\(x)h’(U(x))q(x)  dx  =  A  Cq(x)  dx. 

•'ll  •'(> 

From  (5.6)  and  (5.7)  it  is  clear  that  A  <  0  and  so  U(x)  is  stable. 

When  c(x)h'(u)  >  0.  we  have  two  cases:  (i)  U(x)  contains  no  interior  layer, 
(ii)  L'(.v)  contains  an  interior  layer.  In  case  (i),  for  A  =  A(f)  is  bounded 
uniformly  in  t.  then  we  may  rewrite  (5.5)  as 


(5.8) 


ft/"  ~  (f\U)q)'  +  (c(x)h'(U)  -  A  )q  =  0. 
4(0)  =0.  4'(0)  =  1. 
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which  is  of  the  same  form  as  (4.5).  Thus  it  follow's  from  Lemmas  3.3  and  3.4  that 
q(  1)  >  0.  which  contradicts  (5.5),  <y(l)  =  0.  fn  case  (i)  if  A  =  A(f)  becomes  large 
A  »  1,  then  it  follows  from  integrating  (5.5)  that 

(q(x)E(x))'  =  £  (  V )  -  \e(x)  ['(c(x)h'(U  )  -  \)q(r)dr. 

I  J{) 

E(x)  =  exp |  -  jjf  f’(U )(  t)|  dr: 
whence  (by  vh'(U  )  -  A  <  0  as  A  »  l)we  have 

q{x)>  ih)l'E(T)cJr- 

In  particular.  q(\)  >  0.  again  a  contradiction.  Thus,  in  (i),  A  <  0  and  U(x)  is 
stable. 

Finally  we  consider  the  case  where  c(x)h'(u)  >  0  and  U(x)  contains  an 
interior  layer.  We  want  to  prove  by  contradiction  that  A  >  0.  If  not,  then 

c(x)h'(lt)  -  A  >  0 

and  so  from  adapting  the  proof  of  Lemma  4.5  to  (5.8)  we  have  <y(  1 )  <  0  which 
contradicts  (5.5),  q(l)  -  0.  Thus  A  >  0  and  U(x)  is  unsiahle.  This  proves  the 
theorem. 
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